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ABSTRACT 

It has been established that virial masses for black holes in low-redshift active galaxies can be estimated 
from measurements of the optical continuum strength and the width of the broad H/3 line. Under various 
circumstances, however, both of these quantities can be challenging to measure or can be subject to large 
systematic uncertainties. To mitigate these difficulties, we present a new method for estimating black hole 
masses. From analysis of a new sample of broad-line active galactic nuclei, we find that Ha luminosity 
scales almost linearly with optical continuum luminosity and that a strong correlation exists between 
Ha and H/3 line widths. These two empirical correlations allow us to translate the standard virial mass 
system to a new one based solely on observations of the broad Ha emission line. 

Subject headings: galaxies: active — galaxies: jets — galaxies: nuclei — galaxies: Seyfert — quasars: 
general 



1. INTRODUCTION 

The recent discovery of a tight correlation between black 
hole (BH) mass and the stellar velocity dispersion of the 
bulge of the host galaxy (the Mbh — relation; Geb- 
hardt et al. 2000a; Ferrarese & Merritt 2000) has sparked 
renewed interest in deriving BH masses in active galactic 
nuclei (AGNs). Reverberation mapping (e.g., Blandford & 
McKee 1982) experiments can yield estimates of the radius 
of the broad-line region (BLR) from the lag between the 
variability in the AGN continuum and the corresponding 
variability in the broad permitted lines. With the BLR 
radius (-Rblr) in hand and some measure of the BLR ve- 
locity dispersion (v), one can infer a virial mass for the 
BH, M BH = v 2 R BLR /G (Ho 1999; Wandel, Peterson, & 
Malkan 1999; Kaspi et al. 2000). Remarkably, this simple 
prescription yields BH masses that seem to be accurate to 
a factor of ~3 when compared against the Mbh — c* rela- 
tion (Gebhardt et al. 2000b; Ferrarese et al. 2001; Nelson 
et al. 2004; Onken et al. 2004). The number of AGNs 
with BH masses derived from reverberation mapping is 
currently small (35), but Mbh can be estimated indirectly 
from single-epoch spectra using an empirical correlation 
between -Rblr and optical continuum luminosity that has 
been calibrated using the reverberation-mapped sample: 
-Rblr oc L°^ 00 (Kaspi et al. 2000), where L5100 = ALa at 
A = 5100 A. From a single-epoch spectrum one can mea- 
sure both v and L5100, and thereby infer Mbh- 

While virial estimates of BH masses have been widely 
used in the literature (e.g., Vestergaard 2002; Greene & 
Ho 2004), various systematic uncertainties may affect the 
measurements of both luminosity and velocity. Broad Fe II 
emission in Type 1 AGNs can add ambiguity in the deter- 
mination of the optical continuum luminosity at 5100 A. In 
lower-luminosity AGNs, dilution by host galaxy starlight 
can be so severe that it becomes very difficult, if not 



virtually impossible, to isolate the nonstellar continuum 
unambiguously. As noted by Wu et al. (2004), an ad- 
ditional complication arises in sources, such as blazars, 
whose optical continuum can be subject to contamination 
by synchrotron emission from jets, which might be partly 
beamed. Under these circumstances, the observed £5100 
overestimates the true thermal component that governs 
the BLR size-luminosity relation. 

For nearby AGNs, including the primary reverberation- 
mapped sample (Kaspi et al. 2000), the velocity dispersion 
used in the virial expression is usually measured through 
the width of the broad H/3 line. However, H/3 is at least 
a factor of 3 weaker than Ha, and so from considerations 
of signal-to-noise ratio (S/N) alone Ha, if available, is su- 
perior to H/3. In practice, in some cases Ha may be the 
only detectable broad permitted line in the optical (tra- 
ditionally such objects are known as Seyfert 1.9 galaxies; 
Osterbrock 1981). 

In an effort to minimize uncertainties in virial BH mass 
estimates, and to broaden their applicability to classes of 
AGNs that would otherwise be inaccessible using the con- 
ventional methodology (e.g., blazars, galaxy-dominated 
low-luminosity sources, Type 1.9 Seyferts), we propose to 
translate the virial mass system of Kaspi et al. (2000), 
which makes use of the 5100 A continuum luminosity 
and the broad H/3 line width, to a new system based en- 
tirely on observations of the broad Ha line. We achieve 
this in two steps. First, we revisit the well-known em- 
pirical correlation between optical continuum luminosity 
and Balmer emission-line luminosity (Yee 1980; Shuder 
1981; Ho & Peng 2001), which now can be established 
with much greater precision using large, homogeneous 
samples of AGNs. This allows us to replace L5100 with 
Lh q , which has the advantage of being relatively straight- 
forward to measure even for low-luminosity or galaxy- 
dominated sources (e.g., Ho et al. 1997; Greene & Ho 
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2004). Unlike the optical continuum, the line emission is 
also expected to be little affected by the jet component in 
radio- loud sources or blazars (e.g., Wang, Ho, & Staubert 
2003). This approach is similar to that taken by Wu ct 
al. (2004), who converted £5100 to Lup, but their analysis 
was confined to the original reverberation mapping sam- 
ple of Kaspi et al. (2000), which consists, by necessity, of 
AGNs that are variable on timescales probed by current 
reverberation mapping experiments. This selection effect 
may result in subtle biases, and, in any case, the sample 
suffers from limited statistics. Secondly, we substitute the 
line width of H/3 with that of Ha, using a newly derived 
empirical correlation between the two. 

Throughout we assume the following cosmological 
parameters to calculate distances: Hq = lOO/i = 
71 km s" 1 Mpc" 1 , n m = 0.27, and n A = 0.75 (Spergel 
et al. 2003). 

2. THE SAMPLE 

We utilize the large and homogeneous database pro- 
vided by the Sloan Digital Sky Survey (SDSS; York et al. 
2000), which will eventually obtain imaging and follow- 
up spectroscopy for about one-quarter of the sky. Briefly, 
spectroscopic candidates are selected based on multi-band 
imaging (Strauss et al. 2002) with a drift-scan camera 
(Gunn et al. 1998). A pair of spectrographs is fed by 3"- 
diametcr fibers, covering ~ 3800 - 9200 A with an instru- 
mental resolution of A/AA w 1800 (Gaussian <Ti nst ~ 71 
km s _1 ). Integration times are determined for a minimum 
S/N = 4 at g = 20.2 mag. The spectroscopic pipeline per- 
forms basic image calibrations, as well as spectral extrac- 
tion, sky subtraction, removal of atmospheric absorption 
bands, and wavelength and spectrophotometric calibration 
(Stoughton ct al. 2002). Finally, spectral classification is 
performed using cross-calibration with stellar, emission- 
line, and active galaxy templates. The Third Data Re- 
lease (DR3; Abazajian et al. 2005) of SDSS contains 45,260 
spectroscopically identified AGNs with z < 2.3. Of these, 
~ 3000 have z < 0.35, and hence Ha in the bandpass. 

From this sample we select objects with a high S/N 
(> 25 per pixel) and low galaxy contamination, so that 
we can unambiguously measure £5100- We use the equiva- 
lent width (EW) of the Ca 11 K A3934 line to estimate the 
degree of galaxy contamination; within the SDSS band- 
pass this is the highest-EW stellar absorption feature that 
is not confused with line emission. By artificially dilut- 
ing composite galaxy spectra 1 from Bernardi et al. (2003) 
with varying levels of featureless AGN continuum (which 
for a narrow spectral region is well approximated by a con- 
stant), we find that an EW cutoff of 1.5 A for the Ca 11 K 
line corresponds to a galaxy contribution of < 20%. Thus, 
our primary sample was chosen to satisfy EW(Ca K) < 1.5 
A. In order to increase the dynamic range of the sample at 
the lowest luminosities, we included ~ 30 more galaxies, 
including the 19 objects from Greene & Ho (2004), which 
have a mean S/N w 17 and EW(Ca K) as high as 10 A. 
The total sample consists of 229 objects. 

To investigate the possible influence of jet emission on 
-^5100) w e also consider an auxiliary sample of 59 radio- loud 
objects, selected using the standard criterion for radio- 

lr The composite spectra are constructed from a library of ~ 9000 earl 
rcdshift, luminosity, and effective radius. We degrade the S/N ratio of 



loudness, R > 10, where R is the ratio of the flux densities 
at 6 cm and 4400 A (Kellermann et al. 1989). The ra- 
dio flux densities were obtained from the Faint Images of 
the Radio Sky at Twenty-cm (FIRST; Becker, White, & 
Helfand 1995) survey, which covers all the objects in our 
sample. We corrected the flux densities from 20 cm to 6 cm 
assuming a spectrum /„ cx j/^ - 46 (Ho & Ulvestad 2001). 
The derived R values for the radio-loud sample range from 
10 to 10 4 . 

3. SPECTRAL FITTING 

Below we discuss our method to measure the AGN con- 
tinuum and the broad emission lines. The continuum 
and line measurements are corrected for Galactic extinc- 
tion using the extinction map of Schlcgcl, Finkbcincr, & 
Davis (1998) and the reddening law of Cardelli, Clayton, 
& Mathis (1989). 

3.1. Continuum 

The nonstellar, featureless continuum of Type 1 (broad- 
line) AGNs is challenging to measure at optical wave- 
lengths because it is entangled with a complex emission- 
line spectrum and starlight, depending on the relative 
brightness of the nucleus with respect to the host galaxy. 
Particularly troublesome are the broad, blended Fe II mul- 
tiplets, which form a "pseudo-continuum" that litters most 
of the spectrum, but especially the regions flanking H/3. 
Moreover, the intrinsic shape of the featureless continuum 
itself is not well known. For the objects that have minimal 
starlight contamination, which, by selection, constitute the 
majority of our sample, we decompose the spectrum by si- 
multaneously fitting a two-component model consisting of 
(1) an underlying featureless continuum and (2) an em- 
pirical Fe II template. In practice, the featureless contin- 
uum over the optical region can be approximated well by 
a double power law, with a spectral break at A ~ 5000 A 
(Vandcn Berk et al. 2001; M. Kim & L. C. Ho 2005, in 
preparation). We require that the combined flux of the 
two power-law components at 5600 A, a relatively line- 
free region, be equal to the observed flux at that point. 
Following Boroson & Green (1992), an effective Fe II tem- 
plate can be generated by a simple velocity-broadening and 
scaling of the Fe II spectrum derived from observations of 
the "narrow-line" Type 1 AGN I Zw 1 (kindly provided 
by T. Boroson). In order to avoid contamination from 
strong spectral lines, we only include the following regions 
in our fit: 4170-4260, 4430-4770, 5080-5550, 6050-6200, 
and 6890-7010 A. We do not fit the region below 3685 A, 
since it is not covered by the I Zw 1 Fe II template of 
Boroson & Green (1992); consequently, we do not have to 
model the Balmer continuum (e.g., Wills, Netzer, & Wills 
1985). 

Approximately 30 objects at the lowest luminosities 
have sufficient galaxy contamination to warrant remov- 
ing it. In lieu of more sophisticated methods for starlight 
decomposition (e.g., Greene & Ho 2004, 2005), we find 
that for the present sample the galaxy continuum can be 
modeled adequately using a scaled spectrum of a velocity- 
broadened K giant star (obtained from the SDSS). 

'-type galaxies from SDSS, and span a wide range in velocity dispersion, 
;he composites to match that of our sample. 
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Fig. 1. — Sample spectra for three objects in our sample. The ordinate is in units of 10~ 17 erg s _1 cm~ 2 A -1 . Left: In each panel, 
the top shows the observed spectrum overplotted with the model for the continuum, offset in the ordinate for clarity. The middle shows the 
fitted components (power law, Fe n template, and, if necessary, host galaxy) of the continuum model, with no offset. The bottom shows the 
continuum-subtracted, emission-line spectrum. Center. Fits to the H/3+[Om] AA4959, 5007 region. The data are shown in histogram, the 
fitted broad and narrow components are shown in dotted lines, and the final model (sum of all components) is shown as a solid line. The 
residuals are plotted below. Right: Same as the middle panels, but for the Ho+[N n] AA6548, 6583 complex. 



3.2. Emission Lines 

The primary challenge in measuring accurate fluxes and 
profiles of broad emission lines is to properly deblend con- 
taminating narrow-line components. Here we are only con- 
cerned with the region immediately surrounding Ha and 
H/3. For Ha, as outlined in Greene & Ho (2004; see also 
Ho et al. 1997), we begin by constructing a model for the 
narrow components using the [S II] AA6716, 6731 lines. In 
short, a multi-Gaussian model is fit to the [S II] doublet, 
and then shifted and scaled to fit the Ha+[N II] AA6548, 
6583 complex. The relative positions of the narrow Ha 
and [N II] lines are constrained by their laboratory values, 
and the relative ratio of the two [N II] components is fixed 
to 2.96. We then fit the broad component of Ha with as 
many Gaussian components as needed to provide a sta- 
tistically acceptable fit. This method provides the needed 
flexibility to reproduce highly non-Gaussian, asymmetric 
profiles, but we attach no physical significance to the in- 
dividual components. 

The fitting procedure for H/3 is similar to that outlined 
above; in this case the narrow-line contamination comes 
from the narrow component of H/3 and [O III] AA4959, 
5007. Because of their proximity, we first experimented 
with modeling the narrow H/3 component using [O III] as 
a template, but could not reliably obtain a satisfactory fit. 
The mismatch stems partly from the frequently observed 
asymmetric blue wing of [O III] (Greene & Ho 2005) 2 . The 

2 The blue asymmetric wing of [O m] is presumably absent from the 
asymmetric component (Greene & Ho 2005), fits Ho well. 



profile of [O III] is also expected to be a poor representa- 
tion of the profile of [S II] (and hence narrow Ha and H/3) 
if the narrow-line region is stratified in density (see discus- 
sion in Greene & Ho 2005). Therefore, as with Ha, we use 
the template built from [S II] to model narrow H/3. The 
H/3 centroid is fixed relative to the narrow Ha position 
while its flux is limited by the value appropriate for Case 
B' recombination (Ha = 3.1 H/3; Osterbrock 1989). We 
simultaneously fit [O III] using a two-component Gaussian 
(identical for both lines of the doublet), one for the core 
and another for the blue wing, as described in Greene & 
Ho (2005). Finally, the broad H/3 profile is modeled as 
a multi-component Gaussian, analogous to Ha. While in 
general we find that this procedure provides an acceptable 
fit to the entire region, there were some cases (~ 20) for 
which we could not obtain an acceptable fit for the nar- 
row H/3 using the [S II] model; this may be due to gross 
variations in spectral resolution between the blue and red 
ends of the spectrum. In these cases, in order to fit the 
narrow H/3 component, we retain the flux and centroid 
limits described above but relax the line width constraint, 
and instead impose an (rather arbitrary) upper limit of 
a < 550 km s" 1 . 

3.3. Errors and Uncertainties 

In objects with minimal galaxy contamination, the un- 
certainty in the flux of the optical continuum is dominated 

narrow component of Ha and H/3, since [S n], which generally lacks the 
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Fig. 2. — The correlation between L5100 and (a) Lh„ and (6) I/h/3 f° r the combined broad and narrow components of the line. Filled 
points are our measurements, and open (red) points are from Kaspi et al. (2000). The error bars for the points from Kaspi et al. are smaller 
than the symbol size. The line gives the ordinary least-squares fit to the data. 



by the proper removal of Fe II emission and errors in the 
continuum slope. We use the measured Fe II flux at 5100 
A to estimate the uncertainty contributed by our decom- 
position. To estimate errors resulting from the power-law 
fit, we average the fluxes obtained when the slope of the 
power law is allowed to increase or decrease by one stan- 
dard deviation, as obtained from the fit. The final adopted 
error is the quadrature sum of the Fe II and slope errors, 
although the latter usually dominates. Typical errors are 
~ 0.04 dex. The largest source of uncertainty in mea- 
surements of line flux comes from the placement of the 
continuum level. We estimate our uncertainty by adjust- 
ing the continuum level by one standard deviation above 
and below the nominal best-fit level of the continuum, and 
then recalculating the line flux. Typical errors in the flux 
are ~ 35%, or 0.13 dex. 

Estimating errors in line width is more difficult, since 
uncertainty in the decomposition between broad and nar- 
row components dominates the measurement error. In 
some cases the separation is obvious, but in others (see, 
e.g., the top panel of Fig. 1) there is no clear distinction 
between the two components. We generate an artificial 
spectrum for each object by combining the best-fit model 
emission lines and Gaussian random errors with match- 
ing S/N. An artificial spectrum generated in this manner 
will suffer from the same decomposition ambiguity as the 
original data. We then fit these artificial spectra using 
our fitting procedure, as outlined above. The estimated 
uncertainty is simply taken to be the difference between 
the model and measured FWHM. We find that S/N sig- 
nificantly impacts the reliability of our fits. While typical 
errors are - 7% (or 0.03 dex) and ~ 10% (0.05 dex) for 
FWHM Ha and FWHM H /3 , respectively, these are found to 
increase to as much as 0.1 dex in cases with the lowest 



S/N. 

4. RESULTS 

4.1. Empirical Correlations 

Our data reveal a well-defined correlation between 
Balmer emission-line luminosity and £5100- This is il- 
lustrated in Figure 2, which includes the reverberation- 
mapped sample from Kaspi et al. (2000). Treating L5100 
as the independent variable (under the assumption that 
£5100 traces the continuum that powers the line emission) 
and accounting for errors in both variables (Press et al. 
1992), an ordinary least-squares fit yields 

/j x (1.157±0.005) 

L HQ = (5.25±0.02)xl0- ( lo44 ^ O g ° s _0 CTg S " 

(1) 

and 

/j \ (1.133±0.005) 

L m = (1.425±0.007)xl0 42 ( 1q44 erg S ~\ 

(2) 

The rms scatter in the above luminosity-luminosity rela- 
tions is quite small, being only ~ 0.2 dex in each case. 
The correlation strengths are correspondingly strong, with 
Kendall's t = 0.8, and do not improve when FWHM is 
included as a third parameter using the partial correla- 
tion test of Akritas & Siebert (1996). These relations are 
fit for the combined broad and narrow components. How- 
ever, the best-fit parameters are virtually unchanged when 
only the broad component is considered, since the narrow 
component accounts for only ~ 7% and ~ 10% of the to- 
tal Ha and H/3 flux, respectively Note that the slopes 
of both correlations are formally larger than unity, which 
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can be excluded at a confidence level of > 99% based on 
bootstrap simulations. 
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Fig. 3. — The correlation between Ha and H/3 line widths. Filled 
points arc our measurements, and open points are from Kaspi et al. 
(2000). The solid line gives the ordinary least-squares fit to the 162 
objects with EW(H/3) > 50 A. The dashed line denotes FWHM Hq 
= FWHM H( 3. A typical error bar is shown in the top- left corner. 

We also find a well-defined correlation between the Ha 
and H/3 line width (Fig. 3), which we parameterize 3 by the 
FWHM. In this case, neither FWHM Hq nor FWHM H/ 3 
can be regarded as the independent variable, so we cal- 
culate the ordinary least-squares bisector, again account- 
ing for errors in both parameters (Akritas & Bershady 
1996). Whereas for the line-continuum correlations we 
were forced to include lower-S/N spectra in order to in- 
crease the dynamic range in luminosity, in this case we 
can span the entire range in line width without considering 
the full sample. To minimize errors resulting from profile 
decomposition, we concentrate on spectra with strong H ( 
lines by excluding objects with H/3 EWs less than 50 
For this restricted sample of 162 objects, we find 



Balmer-line luminosity for our auxiliary radio-loud sam- 
ple. If jet contamination is important, we would expect the 
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FWHMh/3 = (1.07±0.07)xl0 3 



FWHM Ha 
10 3 km s- 1 



(1.03±0.03) 



km 

(3) 



The rms scatter around the best-fit line is ~ 0.1 dex. We 
note that when we include the entire sample, or look at 
the radio-loud subset separately, H/3 is still on average 
broader than Ha, but the scatter is larger and the slope 
of the correlation changes slightly. 

4.2. Radio-loud Sources 

If the optical continuum of radio-loud sources is signif- 
icantly boosted by nonthermal emission from a jet, their 
L5100, and hence virial mass, will be systematically over- 
estimated. We can empirically test the magnitude of this 
effect by examining the correlation between L5100 and 

3 Peterson ct al. (2004) advocate the use of the actual line dispersion (second moment of the line profile). We have adopted the FWHM 
here for consistency with the formalism of Kaspi et al. (2000). We also find that the empirical correlation between the line widths of Ha and 
H/3 is significantly tighter using FWHM than the line dispersion. 
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Fig. 4. — The correlation between L5100 and Lh/3 luminosity 
(broad + narrow components) for the radio-loud sources. Filled 
symbols are from our sample; open symbols are from Brotherton 
(1996). Core-dominated sources with core-to-lobe luminosity ratios 
> 1 arc marked as triangles. The line represents our best fit to the 
Lh/3-^5100 relation shown in Figure 26. 

radio-loud objects to lie systematically to the right of the 
best-fit relations shown in Figure 2. However, as shown in 
Figure 4, the radio-loud sources (filled points) do not show 
a significant offset with respect to the main sample. Their 
mean displacement from our best-fit Lu/^-isioo line is only 
+0.02 dex in L51001 and for the strongest radio sources 
(R > 100) the mean displacement is in fact —0.04; neither 
result is statistically significant. Instead of the degree of 
radio-loudness, the more important parameter may be the 
amount of beaming, which can be guaged roughly by the 
degree of core dominance of the radio emission. Unfortu- 
nately, the resolution of FIRST (5" or 10 kpc at z = 0.1) 
is inadequate to distinguish truly core-dominated objects 
from objects with intrinsically compact jets. We thus ad- 
ditionally consider objects from the sample of Brotherton 
1 (1996), which is composed entirely of radio-loud objects 
with tabulated core-to-lobe radio luminosity ratios. For 
the 37 Brotherton objects with measured U-band magni- 
tudes and R values, we calculate £5100 assuming an intrin- 
sic AGN continuum spectrum of /„ oc j/~ - 45 (Vandcn Berk 
et al. 2001). Brotherton's sample is, on average, about an 
order of magnitude more luminous than ours. The entire 
Brotherton sample (open symbols in Fig. 4) is displaced 
by +0.1 dex from the fiducial line of radio-quiet objects, 
increasing to +0.14 dex for the core-dominated systems 
(defined as those with core-to- lobe ratios > 1; open tri- 
angles). We conclude that BH masses based on L5100 
can be systematically overestimated in radio-loud objects, 
but the effect is only serious for relatively luminous, core- 
dominated sources. While the systematic enhancement 
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of the optical continuum relative to Balmer line emission 
in the most radio-loud, core-dominant sources can be at- 
tributed reasonably to jet contamination, we cannot rule 
out the possibility that this effect is due to changes in the 
ionizing spectrum or covering factor of the line-emitting 
gas. 

4.3. A New Formalism for Estimating BH Virial Masses 

The existence of the tight empirical correlations de- 
scribed above (Eqs. 1 and 3) implies that we can trans- 
form the virial mass expression of Kaspi et al. (2000), 
which makes use of L 5W0 and FWHMh/3, into a new for- 
malism that depends solely on the observed properties of 
the broad Ha line (namely Ljja and FWHMh q , derived 
from a single-epoch spectrum), without incurring appre- 
ciable additional uncertainty. 

We begin by rederiving the radius-luminosity relation of 
Kaspi et al. (2000), using the more recent H/3 lags (derived 
from the cross-correlation centroid) and £5100 luminosities 
listed in Peterson et al. (2004), corrected to our adopted 
cosmology. The data points that Peterson et al. deem to 
be uncertain are omitted. We calculate symmetric error 
bars using the mean of the quoted upper and lower uncer- 
tainties, and use the weighted average of all available lag 
measurements. Using an ordinary least-squares fit with 
uncertainties in both parameters, we obtain 



L5100, and for completeness we give the virial mass formula 
based only on the H/3 line: 



M BH = (3.6 ± 0.2) x 10 b 

T N 0.56±0.02 

^H/3 

10 42 erg s- 1 



FWHMh/3 
10 3 km s- 1 



M . (7) 



4.4. Comparison with Kaspi et al. (2005) 

Kaspi et al. (2005) revisit the radius-luminosity relation 
for the reverberation-mapped sample, and we present a 
brief comparison of our respective results. Our best-fit 
-Rblr-^5100 slope, 0.64 ± 0.02, is consistent with the re- 
sult of Kaspi et al., who find a best-fit slope of 0.67 ± 0.07 
for i?BLR derived from H/3 lags alone (line 18 of their Table 
3). The slight discrepancy between our value and theirs 
stems from a combination of the slightly differing cosmolo- 
gies, samples, and fitting techniques employed in the two 
studies. For completeness, we recast our BH mass estima- 
tors based on the radius-luminosity relation of Kaspi et al. 
(2005): 



i? BL R = (20.0±1;1) 



( L 



5100 



V 10 44 erg s- 



0.67±0.07 



It — days. 

(8) 



i?BLR = (30.2 ±1.4) 



/ £5 



100 



\10 44 erg s- 1 



0.64±0.02 



It — days. 
(4) 



Note that the slope of the radius-luminosity relation, 
0.64±0.02, is somewhat shallower than that originally pub- 
lished by Kaspi et al. (2000), 0.700 ± 0.033, although it is 
consistent with the value of 0.67 ± 0.07 recently rederived 
by Kaspi et al. (2005). For v = ufwhm , appropriate 
for random orbits, and using wfwhm = FWHMh/3 for con- 
sistency with Kaspi et al. (2000), the virial formula 4 for 
BH mass is 



Kaspi et al. (2005) consider a radius-luminosity relation 
derived from the average lags of all the Balmer lines, but 
we prefer to focus on the relation based exclusively on 
H/3 lags because the significant difference we find between 
the FWHMs of Ha and H/3 suggests that the lines form, 
at least in part, from kinematically and possibly spatially 
distinct regions. The resulting relations are 



M BH = (1.3 ±0.3) x 10 b 

0.57±0.06 



L 



Ha 



10 42 erg s" 1 



FWHM Ha 
10 3 km s- 1 



2.06±0.06 



M .(9) 



M BH = (4.4 ±0.2) x 10 6 

0.64±0.02 

^5100 



10 44 



erg s" 



FWHMh/3 
10 3 km s- 1 



M Q . (5) 



M BH = (2.4 ±0.3) x 10 6 

0.59±0.06 



10 42 erg s" 1 



FWHMh/3 ' 
10 3 km s- 1 



M Q .(10) 



Substituting Equations 1 and 3 into Equation 5, we obtain 
a virial mass formula that depends on the Ha line alone: 



Aii 



BH 



(2.0t°; 4 ) x 10 6 



U 



10 42 erg s" 



0.55±0.02 



/ FWHM Hq 
\10 3 km s" 1 



2.06±0.06 



M .(6) 



While we argue that Ha is the line of choice whenever 
possible, we recognize that often H/3 is more readily avail- 
able (for instance, due to redshift constraints). In such 
cases, it might still be preferable to use Lh/3 rather than 

4 By comparing reverberation mapping masses with an enlarged sample of bulge velocity dispersions 
normalization factor for the virial formula larger than the one we adopted by a factor of 1.8. 



Kaspi et al. also derive a new -Rblr-£h/3 relation from 
the reverberation-mapped sample. Their slope of 0.69 ± 
0.06 differs somewhat from the 0.56 ± 0.02 found here 
(Eq. 7). While the small difference in adopted cosmology 
plays some role, most of the discrepancy arises from the 
different luminosity ranges covered by our respective sam- 
ples. Our H/3 luminosities range from ~ io 39 - 5 to 10 43 erg 
s _1 , whereas Kaspi et al.'s sample has very few points 
below 10 41 erg s _1 . When use a restricted sample with 
Lh/3 > 10 41 erg s _1 , we obtain a L 5 i o-Lnp relation with 
a shallower slope of 0.97 ±0.04, which would lead to a 
steeper Rblr-Lh/3 relation, consistent with that found by 
Kaspi et al.. A homogeneous sample covering as wide a 

Onkcn ct al. (2004) advocate a 
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luminosity range as possible is required to address the is- 
sue of whether the slope of the -Rblr-^h/3 relation is truly 
luminosity-dependent . 




40 42 



lo S ( L hc / er e s_1 ) 

Fig. 5. — Correlation between Ljj/3 anQl ''Hm f° r the combined 
broad and narrow components of the lines. Filled points are our 
measurements, and open (red) points are from Kaspi et al. (2000); 
radio-loud objects have boxes around them. The line represents 
Ln„ = 3.1 Lh/3i as expected for Case B' recombination. The data 
follow this relation quite closely. A typical error bar is shown in the 
upper-left corner. 

5. DISCUSSION 

The new virial mass formalism based on the Ha emis- 
sion line, which is ultimately calibrated against the mass 
scale of reverberation-mapped AGNs, relies on two strong 
empirical correlations presented in this paper. The first 
correlation, that between optical continuum luminosity 
and Balmer emission-line luminosity, has long been known 
(e.g., Yee & Oke 1978; Yee 1980). Shuder (1981) com- 
piled many of the early measurements, spanning a wide 
range of Ha luminosity from ~ 10 40 to 10 44 erg s _1 . He 
showed that Ln a traces the optical continuum luminos- 
ity with a slope of 1.05 ± 0.03, consistent with unity and 
very similar to the slope found here (1.157 ± 0.005). Simi- 
larly, from an analysis of AGNs ranging from nearby low- 
luminosity Seyferts to quasars, Ho & Peng (2001) found 
that the luminosity of the broad H/3 line, from ~ 10 38 to 
10 44 erg s _1 , correlates strongly with the -B-band contin- 
uum luminosity, roughly of the form Lh/30c L°b 9 - While 
this slope is formally shallower than that found in our 
sample (1.133 ± 0.005), the measurements for the low- 
luminosity sources in the sample of Ho & Peng are also 
considerably more uncertain, making it difficult to judge 
the significance of the apparent discrepancy. Our study, 
based on a large, uniform, and modern data set, provides 
an important confirmation of the earlier results. 

The correlation between the luminosity of the Balmer 
lines and the optical continuum arises naturally if the 
latter traces the low-energy tail of a featureless contin- 
uum that extends into the far- ultraviolet, which powers 
the line emission through photoionization. Indeed, quite 
early on this simple picture provided a natural explana- 



tion for the observed uniformity of the EW of the broad 
H/3 line in AGNs (Searle & Sargent 1968). Recall that the 
slopes of the Balmer line versus continuum luminosity re- 
lations (Eqs. 1 and 2) are formally steeper than unity (1.16 
and 1.13 for Ha and H/3, respectively). This implies that 
the line strength increases mildly with increasing luminos- 
ity [EW(Ha) cx i 5 ioo 016 ; EW(H/3) ex L 5100 13 ], which 
means that the Balmer lines show an inverse Baldwin ef- 
fect (Baldwin 1977). While various studies have found 
a roughly linear relation between optical continuum lu- 
minosity and Lh/3 (e.g., Yee 1980; Shuder 1981; Boro- 
son & Green 1992; Dietrich et al. 2002), Croom et al. 
(2002), consistent with our study, report a positive slope of 
EW(H/3) cx -L5100 ' 19 using composite spectra constructed 
from the 2dF+6dF quasar survey. These apparently con- 
flicting results may be attributable to the different lumi- 
nosity coverage in the different studies. As discussed by 
Croom et al., the luminosity-dependent variation of the 
Balmer line strengths may indicate that the shape of the 
ionizing continuum changes systematically with luminos- 
ity. Another possibility is that the covering factor of the 
line-emitting gas depends on luminosity. 

The second major empirical correlation presented in this 
study links the line widths of the Ha and H/3 lines. While 
a number of studies have investigated the relative differ- 
ences in the profiles of the Ha and H/3 lines (e.g., Os- 
terbrock & Shuder 1982; Shuder 1982, 1984; Crenshaw 
1986; Stirpe 1991; Kollatschny 2003), ours considers by 
far the largest and most homogeneous sample to date. 
Consistent with previous studies, we find that H/3 gen- 
erally has a slightly broader profile than Ha. We measure 
(FWHMh/3 /FWHMhq) = 1.17 with an error in the mean 
of 0.01 and an rms of 0.2. This value is, in fact, identical to 
the average value found by Osterbrock & Shuder (1982), 
but considering our order-of-magnitude increase in sample 
size this agreement must be fortuitous. The relative profile 
of Ha and H/3 depends on both the kinematic and ioniza- 
tion structure of the BLR, since H/3 is emitted preferen- 
tially in regions of higher density and/or higher ionization 
parameter than Ha (e.g. Osterbrock 1989). The tendency 
for the H/3 profile to be broader than Ha is not surprising 
if the density or ionization parameter of the BLR increases 
with decreasing radii. 

Further insights into the physical conditions of the BLR 
may be gleaned from the Balmer decrements that come 
as a by-product of our analysis. The ratio of broad Ha 
to H/3 strengths has been found to vary quite significantly 
from the value of 3.1 predicted by Case B' recombination. 
For instance, while Seyfert 1 galaxies have typical Ha/H/3 
ratios of ^3.5 — 5 (Adams & Weedman 1975; Osterbrock 
1977), in both broad-line radio galaxies and Seyfert 1.9 
galaxies values as high as 7-10 are seen (Osterbrock 1981). 
Interestingly, we find (Ha/H/3) = 3.5 (for the combined 
broad and narrow components; Fig. 5), which is virtually 
identical to the value seen in the SDSS quasar composite 
spectrum of Vanden Berk et al. (2001). [The mean Ha/H/3 
ratio for the radio-loud group is slightly higher (^4), but 
given the smaller sample, this is only a marginally signifi- 
cant difference.] Since this value is only mildly larger than 
the theoretical limit of 3.1, it implies that processes that 
enhance the Balmer decrement, such as collisional excita- 
tion or self-absorption (e.g., Netzer 1975), are generally 
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not that important. We can also place a stringent limit on 
the amount of internal reddening by dust: for an intrinsic 
Balmer decrement of 3.1, the extinction law of Cardelli et 
al. (1989) would allow at most E(B - V) = 0.12 mag. 
Most likely our selection of AGN-dominated sources has 
biased us against objects with large Balmer decrements, 
since the Balmer decrement may steepen in AGNs in the 
low state (e.g., Tran, Osterbrock, & Martel 1992; Korista 
& Goad 2004). 

We will revisit the empirical correlations discussed above 
in a forthcoming study that considers a much larger and 
more statistically complete sample. 

6. SUMMARY 

We present a new method for estimating black hole 
masses from single-epoch optical spectra of AGNs, based 
entirely on the luminosity and line width of the broad Ha 
line. We achieve formal uncertainties comparable to those 
of the previous, widely used method of Kaspi et al. (2000) 
that relics on the AGN optical continuum luminosity and 
the width of the broad H/3 line. The advantage of our 



method is that it enables us to derive robust masses in 
objects where the AGN luminosity may be difficult (e.g., 
jet-dominated systems) or impossible (e.g., low-luminosity 
AGNs with significant galaxy contamination) to measure 
accurately, or in objects where the broad H/3 line may be 
missing altogether (e.g., Seyfert 1.9 galaxies). The robust- 
ness of the Ha-based mass estimator is due to two tight 
empirical correlations we have established: that between 
AGN optical continuum luminosity and Balmer emission- 
line luminosity on the one hand, and that between Hen 
and H/3 line width on the other. We briefly discuss some 
physical ramifications of these empirical correlations. 

We thank T. Boroson for providing us with the I Zw 1 
iron template and an anonymous referee for thoughtful 
comments that improved the manuscript. L. C. H. ac- 
knowledges support by the Carnegie Institution of Wash- 
ington and by NASA grants from the Space Telescope 
Science Institute (operated by AURA, Inc., under NASA 
contract NAS5-26555). We are grateful to the SDSS col- 
laboration for providing the extraordinary database and 
processing tools that made this work possible. 
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